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One of the benchmark analyses to be performed with the first data at the CERN Large Hadron Collider will be 
the measurement of the Z boson transverse momentum spectrum. In this talk, I present a prospective analysis 
for this measurement in the dimuon channel on the ATLAS experiment. The analysis uses simulated datasets 
at a center-of-mass energy of 10 TeV. After summarizing the motivations for the measurement, I discuss the 
Z boson selection criteria, possible physics backgrounds, and background removal techniques with a focus on 
data-driven background determination. I then present my results, and conclude with an outlook toward the 
collision data expected in 2009-10. 



1. Introduction 

The first run of the Large Hadron Collider (LHC) 
at CERN is scheduled to start very soon. A num- 
ber of benchmark analyses are being prepared by AT- 
LAS (A Toroidal LHC Apparatus) and CMS (Com- 
pact Muon Solenoid), the two general-purpose exper- 
iments at the LHC. These analyses will help commis- 
sion the detectors as well as establish known Standard 
Model physics in the new energy regime. 

One of the benchmark analyses is the measurement 
of the transverse momentum spectrum of the Z bo- 
son. The Harvard group is developing an analysis to 
perform this measurement, which is presented in this 
article. In Section [2j we summarize the motivations 
for the measurement. In Section [3l we give the salient 
details of our analysis, including data-driven estima- 
tions of the major backgrounds and background re- 
moval. Section [H presents the results of the analysis. 
We conclude in Section [5j 



2. Motivations for measuring the Z p T 
spectrum 

The transverse momentum spectrum of the Z boson 
is important for a number of reasons, both theoretical 
and practical. The principal among these are briefly 
discussed below. 

• Test of QCD predictions: The high-p^ region of 
the Z pr spectrum is affected by perturbative 
QCD corrections, while non-perturbative correc- 
tions modify the low-px region. As an example, 
Figure [1] shows the Z differential cross-section 
as a function of the Z pr from two Monte Carlo 
programs: Pythia and MC@NLO. Pythia uses 
leading-order matrix elements, while MC@NLO 
uses next-to- leading-order matrix elements. The 
difference in the two distributions is clearly vis- 
ible. 



Recently, a number of Monte Carlo tools have 
become available which combine electroweak 
and QCD corrections. Predictions from the 
state-of-the-art in theory can be tested by com- 
parison with the measured Z pr distribution. 

• Discovery physics: Various physics scenarios be- 
yond the Standard Model predict heavy par- 
ticles that couple to the electroweak sector, 
and therefore can decay to high-pT Z bosons. 
These include technicolor condensates, squarks 
and gluinos from supersymmetric models, right- 
handed quarks etc. If any of these scenarios is 
realized in nature, the decay signature of the 
heavy fields may well show up in the Z pr spec- 
trum. 

• Inferring properties of the Z — > vv decay: The 
invisible Z — > vv decay is a background to many 
processes involving missing transverse energy, 
such as Supersymmetry and extra-dimension 
models where a graviton escapes into the higher- 
dimensional 'bulk'. The Z —> /i/i decay can be 
used to indirectly measure the rate and proper- 
ties of the invisible mode. In order to do this, we 
need to measure the pr spectrum of Z bosons 
in the muon channel. 

• Both the ATLAS and CMS experiments will 
eventually use the total Z production cross- 
section to monitor luminosity on a run-by-run 
basis. Before we can do that, however, we must 
understand the Z pr spectrum very well, espe- 
cially the 1ow-pt region. 



3. Details of the analysis 

Our analysis is based on official ATLAS Monte 
Carlo datasets corresponding to a center-of-mass en- 
ergy of 10 TeV. The datasets were generated with 



2 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



10 = 



> 
O 

■o 



10 1 



i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 
Pythia distribution 
MC@NLO distribution 

ATLAS work in progress 



20 40 60 80 100 120 140 160 180 200 
Zp T [GeV] 



Figure 1: Z differential cross-section vs pr from two dif- 
ferent MC generators: Pythia (black), which uses leading- 
order matrix elements, and MC@NLO (red), which uses 
next-to- leading order matrix elements. 



the Pythia and MC@NLO generators 1 , and were fully 
simulated using the GEANT4 ATLAS detector sim- 
ulation. The sample we use for the signal channel 
corresponds to an integrated luminosity of 40 
which is a realistic estimate of the amount of data 
expected in the first year of LHC running. 

Table U summarizes our Z boson selection criteria. 



Table I The criteria used to select events likely to contain 
Z bosons. The fraction of events passing each cut is shown 
on the right. Errors are statistical only. 



Selection criterion 


Events passing cut (%) 


20 GeV single muon trigger 


69.7 ± 0.4 


At least one /x + and one \i~ 
reconstructed 


49.8 db 0.3 


\n\ < 2.5 for both muons 


47.1 db 0.3 


pr > 20 GeV for both muons 


41.6 ± 0.3 


76 GeV < Mfifj, < 106 GeV 


38.5 db 0.3 



3.1. Background reduction 

Possible physics backgrounds to the signal channel 
include: 



• W — > \iv + jet 

• bb — > fifi + X 

• tt dimuon decays 



• WW, zz. wz 

The first three channels have the largest cross- 
sections, which are the ones we are studying at this 
time. We estimate our backgrounds using Monte 
Carlo samples. However, we are investigating meth- 
ods to extract the W and bb backgrounds from data, 
as we will discuss shortly. To minimize backgrounds, 
we use two isolation criteria defined in a cone of size 
AR = 0.4 around the muon track, where 



AR= a/At? 2 



(1) 



and 77 and (j) are the pseudorapidity and the az- 
imuthal coordinate respectively. The isolation criteria 
are: 

• number of tracks in cone, and 

• the total transverse momentum of tracks, Sp^, 
in cone. 

Figure [2] shows the distribution of these variables 
for the signal channel and the three main background 
channels. We optimize the cuts on these quantities us- 



ing distributions of - 



where5 and B refer respec- 



tively to the signal and the total contribution from all 
three background channels. Figure [3] shows the ~/g=g 
distributions for the two isolation criteria. 

The optimized values of the cuts are given by the 
position of the peaks: 

• number of tracks in cone < 4, 

• YiPt of tracks in cone < 8 GeV. 

Figure [4] shows the distribution of the dimuon in- 
variant mass and the dimuon pr spectra from the sig- 
nal channel and the three background channels. Dis- 
tributions for the quantities are shown before any cuts, 
after the Z boson selection cuts and before the isola- 
tion cuts, and after the isolation cuts, tt dimuon de- 
cays are the largest background contribution after all 
cuts. This is expected, since muons from this source 
are very similar to those from Z decays. Following the 
isolation cuts, the remaining total background con- 
tamination is « 0.2% of the signal, which is small 
enough that we ignore it in the remainder of the anal- 
ysis. 



3.2. Estimation of W 
from data 



[iv background 



1 The PHOTOS tool was interfaced with both Pythia and 
MC@NLO to generate final-state photon radiation. 



In a fraction of the W —> fiu events, the muon from 
the W decay combines with a muon in a jet to mimic 
the Z —> fifi signal. To estimate this background in 
data, we make the following assumption: 
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Figure 2: Number of tracks (top) and Epr of tracks (bottom) in cone of size R = 0.4. Distributions for the muon with 
the larger value of the quantity is shown on the left, and those for the muon with the smaller value on the right. 




Figure 3: Optimization of the muon isolation cuts: (left) number of tracks in cone of size R — 0.4 around the muon 
track, and (right) Y>pt of tracks in cone. On each plot, the arrow indicates the value of the quantity where the peak in 
the distribution occurs. 



The fraction of Z — > events in which three muons 
pass our selection and isolation cuts is equal to the 
fraction of W — » \iv events in which two muons pass 
the same cuts. 

To implement this assumption in a data sample, 
we will first apply the selection and isolation cuts on 



the sample, and then count the number of events in 
which three muons pass all the cuts. Multiplying this 
number by the ratio of the total W cross-section to 
the total Z cross-section will give an estimate of the 
W — > [iv events remaining after the cuts. 

We tested the validity of our assumption on the 
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Figure 4: Distribution of signal and background events in dimuon invariant mass (top) and dimuon pr (bottom) distri- 
butions. Left: before cuts, middle: after Z selection cuts, right: after isolation cuts. 



Monte Carlo samples. In the signal sample, 0.029 ± 
0.008% events (13 events) have a third muon pass- 
ing all the cuts, while 0.027 ± 0.003% events in the 
W — > \iv sample have a second muon passing the cuts. 
These numbers are consistent with each other, sug- 
gesting that our assumption is valid. 

3.3. Estimation of bb background from 
data 

In this case, the idea is to use a sample of bb events 
with non-isolated muons to estimate the number of bb 
events remaining after isolation cuts. In other words, 
we want to use a sample of bb events with non-isolated 
muons as a template for bb events with isolated muons. 
Figure [5] shows the dimuon invariant mass distribution 
in bb events with our isolation cuts and with reversed 
isolation cuts. Reversing the isolation cuts on a data 
sample selects a very pure sample of bb events with 
non-isolated muons. Within statistics, the two invari- 
ant mass distributions have similar shape, indicating 
that the distribution with non-isolated muons can be 
used as a template for bb background after isolation 
cuts. The signal template is a dimuon invariant mass 
distribution from Z — > events after isolation cuts. A 
template fit now gives the fractions in which the sig- 
nal and background templates must be combined to 
obtain the observed distribution. 

We test this method on our Monte Carlo samples. 
The template fit gives 96 ± 1% for signal fraction and 
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Figure 5: Dimuon invariant mass distribution in bb events 
with muon isolation cuts (black) and with reversed isola- 
tion cuts (red), normalized to unit area. 



3.9 ± 0.4% for bb background fraction. The known in- 
puts to our 'data' distribution were 97% signal and 3% 
bb events, so the template fit seems to yield the cor- 
rect numbers. The left-hand plot in Figure [6] shows 
the dimuon invariant mass distributions from 'data' 
and from the template fit result for invariant mass > 
20 GeV; good agreement is observed. The right-hand 
plot in Figure [6] shows the distribution of signal frac- 
tion from the template fit in 1000 pseudo-experiments. 
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Figure 6: (Left) Dimuon invariant mass distribution from data and from the template fit result. (Right) Z fractions from 
the template fit in 1000 pseudo-experiments, in which the content of each data bin was varied using Poisson fluctuation. 



In each experiment, the numbers of signal and bb 
events in each bin of the 'data' distribution were var- 
ied within their Poisson fluctuations, and the template 
fit was done using the same signal and bb templates. 
The width of the Gaussian fit gives an estimate of the 
resolution that can be expected from the template fit. 



4. Results 

The measured Z spectrum must be corrected for 
various detector and reconstruction effects, including: 

• trigger and reconstruction efficiencies, 

• geometric and kinematic acceptance, 

• resolution smearing. 

We derive the correction factors in Z pr bins us- 
ing Monte Carlo truth information (Figure [7j) 2 . Fig- 
ure [8] shows the measured Z pt spectrum from our 
sample before and after the above corrections, as well 
as the true spectrum. Since the correction factors in 
this case were obtained using the true spectrum, the 
good agreement between the corrected measurement 
and the prediction is, of course, expected. 

The differential cross-section as a function of Z pt 
is given by: 



2 With collision data, reconstruction and trigger efficiencies 
will be estimated from data, most likely using a tag-and-probe 
approach. 



dp T * Ap Tz L int €iAi 

where 

Si = signal events in i-th pt bin, 

bi = background events in i-th pt bin, 

Ap^ = width of i-th pt bin, 

Li nt = integrated luminosity, 

Ai = acceptance for i-th bin, 

€i = efficiency for i-th bin, 

Ci = smearing correction for i-th bin. 

Figure [9] shows the Z differential cross-section be- 
fore and after corrections, together with the expected 
distribution from Monte Carlo truth. As before, the 
good agreement between the corrected measurement 
and the prediction is expected. 



5. Conclusion and outlook 

We are developing an analysis to measure the Z 
Pt spectrum in the dimuon channel on the ATLAS 
experiment. During the first LHC run, we expect 
« lOOpfr -1 of delivered integrated luminosity, so that 
we should have a few tens of pb~ x of luminosity on 
tape, assuming realistic data-taking efficiency. Our 
analysis is geared toward that amount of data. We are 
exploring methods to extract the major backgrounds 
to our signal channel with minimal dependence on 
Monte Carlo. 

In the first run, the instantaneous luminosity will be 
in the range 10 29-31 cm" 2 s _1 , so that cavern back- 
ground and event pileup are not likely to be major 
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Figure 7: Corrections to the Z pr spectrum. 




Figure 8: Z pr distribution from Monte Carlo truth 
(black), measured (uncorrected) (blue) and corrected (red). 
The plot corresponds to 40 pb~ x of integrated luminosity. 
Errors are statistical. 



Figure 9: Z differential cross-section vs pr from Monte 
Carlo truth (black), measured (uncorrected) (blue) and 
corrected (red). The plot corresponds to 40 pb -1 of in- 
tegrated luminosity. Errors are statistical. 



issues. The main challenge for this analysis will be 
understanding the detector, including muon trigger 
acceptance and efficiency, offline reconstruction effi- 



systematics associated with all of them. 



